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1 - moment 

The objective* o{ this investigation ate. to mp the variou s thermal 
gradients in tie. c oastal zones o{ France with regard to natural phenomena, 
and man-made thermal eminent* i to study and map the mesoscale thermal 
{eatuA.es in the English Channel f the Bay o{ Biscay and the Month Western 
Mediterranean Sea ; to study and mp the evolution o{ the thermal gradients 
generated by the main estuaries o{ the {rench coastal zones ; and to con- 
tribute to the modelling o{ diurnal heating o{ the sea swi{ace and its 
in{luence on the oceanic sur{ace layers . 

The investigation is conducted by he {ollowings : Dr P.y. PES CHAMPS 
l Principal Investigator.) and p* m. CREPON, Hi J,M. MOMGET and Pro {ess or 
F. MERGER (Co-Investigators J . 

Appendix. A give related organizations and addresses. 

This progress report is the last one be{ore {Inal report o{ the 
investigations , Results have been emphasized, while methods and problems 
have not been discussed. 

2 - RESULTS 

2.1. Residual { low through the Dover Strait 

Time sequence o{ HCMM scenes allowed us to outline the in{luence o{ 
meteorological conditions on the residual current which {lorn to the N.E. , 
^a.oih the British Channel into the Month Sea, through the Dover Strait.. 

S. W. winds enhance this residual {low, and, as a result, the thermal 
e{{luent o{ the Rhine River is kept to the Dutch coast in a very narrow 
coastal band, M. E. winds contrary the residual {low which is slown down 
and deviated to the English coast : then the Rhine thermal e{{Zuent pro- 
pagates o {{shore at a distance o{ up to 25 nautical miles. A close corre- 
lation exists between wind speed direction and the o{{shore width o{ the 
e{{luent. 


2.2. Upwelling at the conttnentol. skelb bleak in the Bom ob tiscau 


HCMM data, conbim the. existence ob a pemancnt upwelling at the con- 
tinental shelb bleak in the Bay ob Biscay. The upwelling it outlined by 
the appeaiance ob c old watei in sumeitime. This was pieviously menti.onned 
biom MRS data. Tim HCMM scenes, a moie emplete description and intelpie- 
tation ob the upwetlirig hat been obtained. 

( I) The, upwelling it piobaliy peimanent, but it enhanced by up- 
welled coldeiwateiinswmeitimewhenaseasonaiiheAmocline is { oimed . On one 
occation, januaiy 16, 1 979, waimei wa tei appealed in wintentime at the 
shelb bleak (HCMM 6 cent n° 265 - 1090 ) .* this watei it piobally a "meditei- 
lanean 1 ' watei, women and tatted, blowing out ob the mediteiianean tea, 
biom the Gibialtai Stiait, into the Atlantic, at a depUi ob several hun- 
ched meters. 


(2) Upwelling it enhanced abtei spring tides, which suggests 
that the basic mechanism joi the upwelling it a tidal one. On two occasions 
abtei spring tides, august 25 and September 21, 1978, HCMM scenes 
(IP n fl A-A0121 - 15260 and A-A 0148 - 13320) show very similar patterns ob 
cold watei at the shelb bleak, with a maxunum intensity between 48H-8E and 
46*30N-5E where the tidal currents aie at a maximum. 

2.3. lietoscale variability ob the SST bield . 


Using VH RR and HCMR inbiaied digital data, a statistical two-dimensional 
analysis ob the mesoscale variability ob the SST bield has been peiboime.d in 
oiden to chaiacteiistize the landom piopeities ob thit bield. The powei law 
exponent, a, ob the spatial variance density spectium, E ( fe) ~ k~ n fe is wave- 
number) , is deduced biom the computation ob the stiuctuie banction ob the 
SST. The study was birst stalled on MHRR/N0AA-5 in the lange ob scales 40« 

100 km. HCMR data allowed us to extend the study down to a scale ob 5 am. In 
the lange og' scales 3-100 km, n was gound to vary biom 1,5 to 2.3, with a mean 
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value oh 1.8, over a study otf 11 WRR and. 9 HCMM <ic enes. These values oh n 
are oh the order oh the. predicted values by the two-dimensional turbulence 
theories . However a discrepancy exists and we need burthen advanced theories 
to explain this experimental determination oh the mesoseale SST variability. 

The. inability oh the spectral analysis In the range oh 4 cole* 3-30 km 
was t made possible by the only low noise level oh the HOAR data. A detailed 
manuscript is given as AppencUce B. 

2.4 . VuirnaR. heat ing 

Vaytime HCMR data occasionnaly exhibit warmer tea surhace areas which, 
extend over 10 to 100 km, The wanning .is oh several 9 C and it eas-iiy detected 
on' photographic product* because the warmer area* have usually smooth boon - 
dories and cannot be conhus ed with the sharper oceanic thermal boundaries . 

These warmer areas are Interpretated as a large diurnal heating oh 
-the upper surhace layer under low wind speed conditions , Evidence oh that Is 
supported by several arguments. 

|J) Meteorological observations and analysis show that warmer areas are asso- 
ciated with low wind Speed conditions - l.e. antlcyelonlc conditions or coastal 
breeze eh heats. 

12) Glitter - l.e . direct solar radiation rehlected by the wavy sea surhace 
towards the sensor - has been used to derive an equivalent wind speed hrom 
the HCMR visible channel, where h^rslble [observation must be dose to the 
specular rehlection otf a hint sea) ■Warmer areas are always associated with 
changes in the glitter patterns and decreasing wind speeds. 

(3) Warmer areas disappear on consecetuve nlghtime HCMR data. 

Under these low wind speed conditions, ZunJbulence Induced in the sur- 
hace layer by the wind stress is strongly reduced, and most oh the solar ra- 
diation absorbed Is stored', without downwards propagation. Theoretical simu- 
~ lotions using a radiative and heat irons her model have been perhonmed and 
predict large heating rates in the upper meter, and a maximum heating oh 
several e C in the upper layer which Is conhitmed by a hew in-sltu measurements. 
Large heating only occurs In a hew tens oh cw and is very rapidly destroyed 
by the nlghtime cooling . 


HCM R data, a (lowed us to discover that a diaanal heating oh moat than 1* C 
could a&itct laage onto*. frequencies oh oacuaenco. oat xelativity high in the 
Mutton M editexoanean Sta wheat moat than 10 $ oh marine suahaae ant abet- 
ted one day on an other, white a laagt diunnaJL heating is very unlikely in 
the Hoath Sea {only one Acene). In such strongly abetted anew, daytime satel- 
lite data could consequently give meaningless SST fields, and observations 
should be xestxicted to nightime, ox eaaty in the moaning when the suahace 
layea is the nhst homogeneous, A detailed manuscript is given as Appendice C. 

3 * COUCLUSim 

Vuaing the aepoating peaiod, HCMM photogxaphic poo ducts paoved to be 
veay useful : 

(1) to inteapxet the infituence orf wind diaection on the mean xesidual Alow 
thaough the V oven Stoait, 

(2) to understand the upwelling occuoing at the continental shelh baeak in 
the Bay oh Biscay, and its oetation with tidal cuxxents. 

A multitempooal analysis HCMM digital paoducts is on pxogxess 

(?) to obtain a mean value oh diuonal heatings obsexved in the Western 
Meditexxanean Sea, duoing surnex months, 

it) to obtain a quantitative assessment oh the intensity oh the shelh 
baeak upwelling in the Bay oh Biscay as function oh tidal conditions. 
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ABSTRACT 


Satellite Infrared data have been used to inveatigate th<? mesoscale 

j 

variability of the S ST (Sea surface Temperature) field. A statistical two- 

i 

dimensional analysis of the SST field has been performed by means of the I 
structure function. Results give the equivalent power law exponent, n, of J 

r'A 

-n ^ 

the spatial variance density spectrum, E(k) *v k . n was found to vary from 
1.5 to 2.3 with a mean value of 1,8 in the range of scales 3-100 km# in , 
agreement with-. previous one-dimensional analysis from shipbome and airborne 
measurements . These observed values of n are discussed and compared to the 
values predicted by turbulence theories. 


1 ( introduction 


The present capability of satellite infrared radiometers permits tha 

determined; on of the mesoscale SST (Saa Surface Tempera tur a) field on an 

operational basis thanks to their improved radiometric performances, which 

2 

are typically a few tenths of (*C) for a nadir resolution of 1 km . This 
gives a potential tooil for a systematic investigation of mesoscale thermal 
features such as thermal fronts and gradients which have already been detected 
and studied by means of infrared imageries or derived SST maps. Besides these 
observable features, a part of the SST field oust be considered as random 
and containing some other information which can only be retrieved by a sta- 
tistical analysis - e.g. the spectral density g£ variance, 

Attempt* to compute the spatial spectrum of the SST have been 
previously made by Me Laish (1970) , Saunders (1972 a) and Holladay and 
O'Brien (1975), from airborne Infrared measurements along the aircraft 
track. Examples of mesoscale spectra have also been given from shlpbcme 
measurements (Voorhis and Perkins , 1966,. Fieux et al. , 1978), but more fre- 
quently for time series than for spatial variations. On the other hand, the 
satellite observations give the unique opportunity of investigating the me- 
soscale variability of the SST field, in the two-dimensions, down to scales 
of 1 km, at any given time, with a frequency which is limited only by the cloud 
cover. In this study, we intend to demonstrate the feasability of using sa- 
tellite data to obtain statistical parameters of the mesoscale SST field. 

2. Statistical analysis of the SST field 

Studies of the variability of the temperature (or any scalar) field 


usually tasks an extensive use of spectral methods * l.e. the computation of 
the spectrum of the density of the scalar variance by means of Fourier trans- 
form or autocorrelation function , to obtain a typical power law which 
characterizes the variability of the temperature field and which cm be 

referred to turbulence theories. In the present study, the structure function 
has been preferably used to determinate more accurately the power law exponent 

in the presence of the large noise level of satellite infrared data. 

* 

2.1 - Structure function 


The SST field is considered as an isotropic random process with homo- 
geneous increments (locally homogeneous) for which the structure function 
can be cuv-iutea as : 

D (h) 4 E (T(x+h) - T (x) } 2 (1) 

IT * 

where T(x) « temperature at x, 

■> 1 

h » scale , k « h ■ wavenumber , 

E •* average operator . 

The main advantage of the structure function, D (h) when compared to the 
spectrum of the variance density, E(k) , or the autocorrelation function, 

B (h) , is that the experimental determination of the structure function is 
more accurate and much less affected by random variations because only incre- 
ments are taken into account (Panchev, 1971) . An example is given in Fig, 1 
where both E T (k) and D^h) have been computed and are shown for the same 
sample of the SST field, measured by the AVHRR (Advanced Very High Resolution 
Radiometer) experiment on board the TXROS-N satellite. This example shows’ * A * v 


clearly that the structure function is more regular than the spectrum, allowing 
an easier determination of the characteristic parameters - e.g. the 
power law exponent given by the slope when using logarithmic coordinates. 

2..P -Interpretation of the structure function 

The structure function D(h) statistically represents the influence of 
a point upon the h- distant points. For an homogeneous and isotropic random 
process > D(h) and B(h) are linked by the following relationship : 


D(h) ■ B (o) - B(h) 


( 2 ) 


As B(h) and E(k) are the Fourier transforms of each other, D(h) may thus be 
related to E(k) (Panchev, 1971) : 


D(h) ■ 2 


( 1-J Q (kh)) E(k) dk (2-dimension, field) 


(3a) 


D(h) « 4 


(1 - E (k) dk (3-dimension field) 


kh 


(3b) 


where J (kh) is the zero order Bessel function of the first kind, 
o 


In geophysics , the spectrum is usually expressed in the following way ; 


E(k) « B k 


-n 


(4) 


and then, from (3a) , (3b) , the structure function may be written as 


D(h) = A h 


(5) 


£ - 6 - 


where A and B are constants , 

n and p are exponents, which may be deduced from each other by ; 

n * p + 1 . (6) 

so that the exponent, n, of the spectrum can be altarnately determined from the 
structure function, using (6) , as far as the field under study is homogeneous. 

m 

.Two kinds of error may affect the determination of the SST field 

* 

obtained from satellite ? instrumental data noise, atmospheric effect. 


Although the structure function has the advantage of being much .more 
regular them the spectrum, the study of the structure function and of its shape 
is generally limited by the noise level at the smallest scales. This effect 
is illustrated in Fig. lb, where the observed slope giving the power law ex- 
ponent of the structure function decreases from about 1 at larger scales, to 
O at smaller scales. This is due to the fact that the structure function of 
the instrumental noise adds to the SST one . As far as this noise is white, 
its structure function is a constant (p«o) and its addition restricts the sta- 
tistical analysis at smaller scales. This effect can be reduced by spatial 
smoothing with a corresponding degradation of the ground resolution. 

Smoothing also introduces a bias in the determination of the structure 
function. If D^th) is the structure function of the Smoothed field, and Q is 
the convolution square of the smoothing function F, it may be shown (Matheron, 
1970) that : 


D (h) - D » Q - A 


( 7 ) 


where §, neons convolution end A is e constant,: 


A * 


f 

+® 

D (U) Q (u) du 


(8) 


In the particular case where F is the spatial average in/'square and where the 
structure function D(h) h^,with 0.5<p< 1.5, the influence of smoothing on 
the amplitude of the structure function o (h) increases with p^ but decreases 
rapidly when h Increases ,and is less than 10 % when h is equal to 5 times the 

9 

dimension of the smoothing square. The influence of spatial smoothing was thus 
neglected in the present study. 

As far as the variations of the atmosphere can be neglected within the 
mesoscale oceanic range, the observed satellite variations of the SST field 
are reduced by the atmospheric infrared transmittance, t (Deschamp6 and 
Phulpin , 1980) : 

T(x+h) - T(x) ■ t (T g (x+h) - T g (x)) ( 9 ) 1 


where T g is the actual SST, 

T is the measured SST from space. 

ftien : 


D (h) 

TO 


C lh > 

TO 


( 10 ) 


where is the actual structure function, x depends on the atmospheric water 

vapor content and ranges between typical values of 0.9 to 0.3 for the 10.5 - 
12.5 pm channel mostly used on satellites. This atmospheric effect will affect 

the determination of the amplitude of the structure function, but not the 
determination of the power law exponent, p. Because the atmospheric transmit- 
tance cannot be accurately determined over the oceans, only one parameter of 
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the structure function can be determined from satellite, and this is the 
power law exponent, p, obtained from the slope of the curve in a log~log plot. 

The hypothesis of the homogeneity of the random field must be verified, 

Otherwise erroneous determinations of the exponent could be obtained. For 

-2 

example, a frontal rone would have a spectrum £ (k) n, k and a structure 

_2 

function D (h) 'vh, while a non-removed trend would also give E_(k) n, k , 

TT 1 * 

2 

but D TT (h) 'v k . Because these exponents are close to the values physically 
expected, it is necessary to check earefully the homogeneity of the SST field 
and to remove the existing trend when necessary. 

3. Results. 

The results of two independent but complementary studies are hereby 
presented. The first one deals with data obtained from the VHRK (Very High 
Resolution Radiometer) on board NOAA-5 and the study was limited to the range 
of scales 40-100 km because of the large level of instrumental noise. The 
improved radiometric performances of the HCMM (Heat Capacity Mapping Mission) 
data, - i.e. a nadir resolution of 0.5 km and NEDT ■ 0.3° K(see Table 1) - 
allowed us to extend the study down to scales of 3 km. 

Cloudfree satellite data were selected in homogeneous study areas, 
Northeastern Atlantic . ocean and Mediterranean Sea. Locations are shown in 
Fig. 2 and dates are given in Table 2. At each one of these locations, the 
unidimensional structure functions were computed in four directions, (across 
the satellite track - i.e. approximatly east to west), ir/4 , ir/2 (along the 
satellite track) and 3ir/4. More details on the processing of the data may be 


found in Frouin ( i960) and Hald ( 1980) . Examples of the computed structure 
functions are given in Fig. 3 for VHRR/NOAA-5 and in Fig, 4 for HCMM. The 
results generally show that the SST field is not exactly isotropic. Nether* 
theless , the structure functions , if not equal, are roughly parallel on a 
log- log plot, so that the anisotropy is confined in the amplitude, A(0) : 

D (0,h) - A(0) h p (11) 

TT 

but the slope p remains very isotropic. 

Values of p from 0.5 to 1.3 have been observed in this study with an 
estimated accuracy of about o.l. Using VHRR/NOAA-5 data, 44 estimations of p 
were made in the range of’ scales 40-100 -km, and 37 estimations in the range 
of scales 3-30 km with HCMM data. The corresponding histograms of the obser- 
ved p are given in Fig . 5a and Fig. 5b. The most frequent values are 0. 9-1.0 
and the mean values are 0.8 (3-30 km) and 0.9 (40-100 km) with a standard 
deviation of about 0.2. About 90 % of the observed values are distributed 
between 0.5 and 1.1. The results correspond to a mean value of the power law 
exponent of the spectrum, n, of 1.8 in the wavenumber range 0.01-0.3 km . 

-2 ui 2 

The amplitude of the structure functions varied from 10 to 10 (*C) 

at h « 40 km. Even after spatial smoothing ,it was noted that the noise level 

had a slight tendency of reducing the estimated of values of p because the 

structure function of the noise is a constant (p«o) . This is particularly 

effective when the noise level (5 . lo”"^ (°C.) ^ for the HCMM data, 3. ,10 ^ ( °C.)^ 

for the VHRR/NOAA-5 after smoothings) is of the same order as . the structure 

function (see Fig. 1) . Whenever possible, the estimates of p were corrected 

for this effect, but the effect could partly explain the lowest values of p. 


On the other hand, a mean horizontal thermal gradient would give 

2 

D(h) a* h . The areas studied were carefully selected to avoid the existence 
of such thermal gradients which would increase the estimate of p towards 
larger values ; but here again some influence on the data could remain. Both 
these two effects, noise level and horizontal thermal gradients, could par- 
tly but not totally explain the spread of the results around the mean value , 
between 0.5 and 1.3, which remains significant. There is no evidence of cor- 
relation between the estimates of p and the corresponding geographical loca-* 

tions or seasons,' but nevertheless, we would guess that it is probably necessary 

* 

to involve physical processes in the explanation of the observed p values. 

4, Discussion 

Using Eq. 6 and the result from thi s structure function analysis, we 
obtain 1.5<n < 2.3 for the power exponent of the spectrum. This agrees 
fairly well with the previous results reported by several authors either 
from shipborne measurements (Fieux et al, 1978) , 

or from airoorne infrared measurements (Saunders , 1972a) , for the one-dimen- 
sional temperature spectra (see Table 3 ). Holladay and O'Brien (1975) also 
made an attempt to reconstruct 'idle two dimensional SST field from the tracks 

* 

of the aircraft survey and found n * 3. for the isotropic part of the two-di- 
mensional spectrum, a value which is probably overestimated because of the * 
smoothing of high wavenumbers produced by the SST mapping procedure . 

it would be interesting to relate the computed values of n to those 
given by turbulence theories in geophysics . All these theories assume the 
existence of an inertial range, i.e. the considered^Sf^ e far from the energy 
sink and source scales. It is not evident that the range of scales 3-100 km 
in the ocean is an inertial one. The upper limit of the dissipation scale 


is of the order of 500 m and can be related to the wavelength of surface 
and internal waves via breaking processes. This scale is about one decade 
smaller than the lower limit of the studied range and we consider that there 
is no interaction between these two scales. The scales of imput of kinetic 
energy and of temperature variance remain puzzling. Input of kinetic energy 
related to the wind i§ found at scales of the order of the wind waves 
(100 m) and at scales of the meteorological systems (1000 km or more) . Energy 
inflow due to thermodynamic forcing is found at even larger scales. All 
these scales are one or two order of magnitude smaller or greater than those 
studied. At some location , interior processes such as baroclinic eddies or 
baroclinic instability may also play an important role in converting energy 
through non-linear mechanisms. The scales of these phenomena are of the same or- 
der as the internal radios of deformation or two to six times greater# depen- 
ding on initial conditions. This radius is^fbout 50 km in the open ocean and 
10 km in the Mediterranean- sea. If these physical processes are of importance 
in the area studied# the range 3-100 km is not an inertial one. In fact# we 
cannot precisely determine this from our observations : by looking at 

Fig. 3 and 4# one can notice that the structure functions do not exhibit any 
peak characterizing a very^§§ffe t rn the range we deal with, but this may 
only mean that the energy inputs are outside the studied range. 

In the range of scales 3-100 km, horizontal scales are larger compared 
to vertical ones and the observed variability may be considered as being a qua- 
si two-dimensional process. Therefore the observations can be related to the 
n- values predicted by the theories of bidimensional turbulence (Kraicftiian , 1971) 
and of geos trophic turbulence in the atmosphere (Chaxney , 1971). These theories 
take into account either the conservation of energy and the conservation of 
enstrophy (half of the mean square of the vorticity) in the case of Kraichnan's 
theory or the pseudo-potential enstrophy (Chamey) . Both these theories agree 
when predicting the power law of the kinetic energy spectrum : (k) ^ k~ 3 . 


But the relations between current and temperature are not Obvious and the 
different mechanisms involved lead to drastically different theoretical power 
laws for the temperature variance spectrum. In Krai.chn.an 1 s theory, consi- 
dering that temperature is a passive contaminant implies that (k) only 
depends upon k and upon the dissipation rates of enstrophy and temperature 

variance. Then. from a dimensional analysis, £ (k) must follow a k”* power 
' T 

law. Chamey made use of the perfect ga^Law and of the hydrostatic relation to 

compute a relation between the temperature and the stream function and .he 

-3 

found the same law for E (k) as for E-Jk) - i.e E (k) ^ EJk) ^ k . Further- 

T ' K t K 
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more, assuming also E^(k) k , Saunders (1972b) exhibited a temperature 

-5 

variance spectrum E^(k) ^ k , by the use: of the thermal wind equation. These 
examples show that the results may be very different according to various 
authors. In this study, the mean observed value of 1.8 for n is far from the 
assessment (n«5) made by Saunters but falls between the Kraichnan and Chamey 
predictions (nsl and 3) . This discrepancy may be due to the fact that the 
theories hypothesises have not been respected,^ particular the hypothesis that 
the range studied is not an inertial one. 

Some three- dimensioned, theories of turbulence (Kolmogorov, 1941, 
Bolgiano, 1962) or space-time variability theories of internal waves (Garrett 
and Munk, 1972, 1975) report values of n close to those found in our 
study (respectively 1.7, 1.4 and 2) , but their physical basic hypothesis can 
hardly be extended in the mesoscale range. 

We may also notice that several studies of atmospheric temperature 
fields mention values of n in agreement with our study at similar range of 
scales (100-1000 km). One may refer to the reviews by Gage (1979) and 
Panchev (1971). some of these results are obtained by using spectral analysis 


on time-series data and equivalent wavenumbers are computed by using Taylor's 
relation. As the validity of this relation is dubious for such scales, these 
time-series results must be considered carefully. But as for the oceanographic 
observations, there is no atmospheric theory to explain the observed results. 

In summary, the power law exponent n of the spectral temperature variance 
observed in the range of scales 3-100 km is nearly 2. a large discrepancy exists 
with the predicted values from the 2 - dimensional turbulence theories which 
are widely spread around this value, and we need further advanced theories 
to explain the experimental determination of the mesoscale SST variability. 

5, Conclusion 

In this study, it has been proved feasible to estimate the random pro- 
perties of the SST field in the mesoscale range 3-100 km from satellite in- 
frared data. Compared to previous 1-dimension analysis from shipborne and 
airborne observations, the use of satellite data allowed us to perform a 2 - 
dimensional analysis. Using the structure function, the power law exponent, 
n, of the spectrum of the variance density of the SST field can be retrieved 
within a good accuracy (+ 0.1). A mean value of 1.8 and a standard deviation 
of 0.2 have been found in the range 3-100 km, and extreme values of 1.5 and 
2.3 have been observed. 

The results give rise to several questions : (.j) Is the range 3-100 km 
an inertial one ? (ii) If yes, is there any turbulence theory to explain the 

•if’--- 

spectrum power law observed ? (iii) If not, at which scales are the inputs of 
energy and to which processes are they related ? At the present time,, further 
investigations, both theoretical and experimental, are needed to interpretate 


the physical mechanisms and paramatars involved in the mesoscale variability 
of the SST field, 
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Table 2, Summary of the different areas studied 

Area 

i 

Date 

Location 

Experiment 


Eastern Mediterranean Sea 

19 Mar*, 1978 

33°00 , N-28*00*E 

VHRR 


I* 

05 May ,1978 

34*00 , N-15 f Q0'E 



H 

08 May ,1978 

33*00*N-29*00'E 

■M 


II 

14 May ,1978 

33*30'N-28 # 30 , E 

H 


II 

17 May ,1978 

33 8 30*N-26 # 00'E 

It 

* 

Western Mediterranean Sea 

29 Sep, ,1977 

41 8 00*N«04 8 00'E 

•I 


II 

29 May ,1978 

39°Q5'N-07 8 15'E 

HCMm 

# 

II 

29 May ,1978 

40 8 05'N-06 8 55'E 

II 


If 

11 Jul.,1978 

38°55'N-04°50 , E 

If 


II 

11 Jul.,1978 

41 8 55 , N-06 8 55'E 

II 


II 

26 Jul.,1978 

39 8 20'N-06°15'E 

II 


If 

28 Jul.,1978 

38°15'N-03°45'E 

II 


II 

28 Jul.,1978 

38°35'N-05 8 05'E 

II 


II 

28 Jul.,1978 

37 8 40*N-07 8 25'E 

II 


ft 

14 Aug., 1978 

38°30'N-03°00'E 

VHRR 


II 

14 Sep, , 1978 

40 8 25 ' N-06 0 30 ' E 

ncm 


li 

1. Sep., 1978 

40 8 35'N-11°55'E 

I* 


II 

14 Sep., 1978 

41 8 40 , N-06 8 45'E 

ii 


Northeastern Atlantic Ocean 

11 Sep,, 1977 

46 8 00'N-06 8 30'W 

VHRR 

s 

it 

14 Sep., 1977 

45°00‘N-07 8 00’W 

ii 


ii 

06 Jan., 1978 

46 8 30 , N-09 8 00'W 

ii 


n 

10 May ,1978 

46 8 00'N-08 8 00'W 

if 


it 

11 May ,1978 

45°15'N-04 8 40'W 

HCMM 


ii 

11 May ,1978 

38°35 'N-ll °45 ' W 

II 


ii 

18 Jun. ,1978: 

46°00'N-08°35'W 

tf 




Table 3. Summary of observed mesoscale SST variability. 




CAPTIONS 


Figure 1 - Comparison between the density of temperature variance (k) (a) 

and the structure function D. (h) (b) , computed from AVHRR data, 

IT 

July 17, 1979, over the Bay of Biscay (45° 30* N - 4° 30* W) . 

The dashed line indicates the radiometer noise level. 

Figure 2 - Geographical locations of the different study areas for HCMM 
data (e) , and VHRR data (•) . 

A 

Figure 3 - Example of structure functions computed from VHRR data. 

Figure 4 - Example of structure functions computed from HCMM data. 

Figure 5 - Histograms of the observed values of the power law exponent p of 

the structure function in the range of scales 40 - 100 km (a) 
and in the range of scales 3-30 km (b) . 






data (<*) , and VHRR data (•) » 
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ABSTRACT 


Day-night surface temperature differences have been measured in the 
infrared (10.5 - 12.5 pm channel) by the HCMM satellite experiment, which 
show large diurnal heating (several 9 C) of the upper layer of the ocean, 
very frequently during summer months in the Mediterranean Sea, when the wind 
speed is low. When observed in the 0.5 - 1.1 pm channel, glitter reflectance - 
i.e. direct solar radiation specularly reflected towards the sensor - correlates 
whith diurnal heating. Glitter reflectance has been modelized to retrieve an 
equivalent wind speed , and observed diurnal heatings, AT, rapidly decrease with * 
the wind speed, U, from a maximum value of about 5° C. An empirical law is 

-3 

given : A T (o K ) - 3.5.10 2(w , m ~ 2 ) / (0.7 + s _ij) where Q is the Irradiance 

at sea level. A mean diurnal heating of nearly 1° C is calculated for the 
marine coastal areas of the south France, in summer time. During this period, 
satellite observations should be restricted to night and early morning times, 
and to the only high wind speed (U > 5 m.s *) . u 
afternoon 


at noon and during the 


I - INTRODUCTION 


A daily variation of the temperature in the surface layer of the oceans 

is known to be produced by the diurnal heating of the absorbed solar radiation. 

The amplitude of the daily temperature is usually small because of the large 

turbulent mixing which usually prevails over the molecular thermal diffusivity. 

_2 

A solar irradiance of 1000 w.m when absorbed in a mixed layer of 10 m would 
only give a heating rate of o.l a C per hour, and a daily variation of less 
than 0.5 9 C. While if the turbulent mixing is reduced and the mixed layer 
thickness is restricted to less than 1 m, a heating rate of 1° C per hour 
may be expected and daily variations of several °C should be observed. With 
the exception of very shallow waters, large diurnal heatings in open oceans 
thus correspond to the case of lower wind speeds as far as turbulence in the 
upper surface layer is mostly locally in".duced by the surface wind stress. 

Prom a theoretical simulation of radiative and heat transfer 
in the upper ocean layer, BASSE (1971) has predicted the deviation of the sea 
surface temperature (SST) T q from the bulk temperature T 1Q taken at 10 meter- 
depth should vary as : 


T - 

o 


‘10 


- C 2 Q 0 


-1 


(1) 


-3 

where Q is the solar irradiance, U, the wind speed, and C. = 3.5 10 when 

-2 _i 

Q is expressed in W.m” , U in m.s” . According to , Eq . ilj is only 

valid for U >, 2m. s * , but the evidence that the SST diurnal variations incre- 
ases when U decreases is supported by several observations : ROMER (1969), 
STOMMEL et al (1970) where large diurnal variations of more than I s C are 
occasionaly found at very low wind speeds - i.e. for U < 2m. s”*. These obser- 


vations art nevertheless rtttricttd to a single location and limited time 
occasions . 

Satellite infrared radiometers offer the opportunity to investigate 
more systematically such large diurnal variations of the SST, The first sa- 
tellite experiment to provide adequate capability for this purpose was the 
HCMR (Heat Capacity Mapping Radiometer) experiment launched in late 
April 78 with an improved temperature resolution (0.3° C) and a nearly noon 
overpass. Results from this experiment are hereby reported in order (i) to 
investigate large diurnal SST variations at low wind speeds (ii) to give an 
assessment of the relation frequency of such an event and its impact on the 
determination of the SST field in such area -W* Mediterranean Sea where the 
occurence of diurnal heating is rather large. 

II . OBSERVATIONS OF DIURNAL HEATING FROM HCMR SATELLITE DATA 

II- 1 - The HCMR experiment 

The basic objectives of the HCMR experiment are to measure diurnal 
variations of the earth surface temperature for applications to earth re- 
sources (geology, hydrology...). For this purpose, the satellite is sun-syn- 
Ch ' * : 

chronous and orbit was chosen to cross the equator at about 2 a.m and 2 p.m 

local time so that surface temperature data are obtained close to the minimum 
and the maximum of the diurnal variation. Satellite altitude is 620 km, and 
orbit inclinaison is 96.87° . The HCMR consists of a two-channel scanning 

radiometer, with a 0.5 - 1.1 ym spectral bandwith in the visible and 10.5 - 
12.5 ym in the thermal infrared. Similar channels have been used on previous 
meteorological satellites, but the interests of the HCMR experiment are (i) a 


large improvment of the radiometric performance* in the thermal infrared 
channel for which the temperature resolution is 0.3* C and the nadir ground 
resolution is 500 m as compared to respectively 0.7 0 c and 1 km for the pre- 
vious VHRR/NOAA satellite, (ii) the facility offered to the user to obtain 
differential surface temperature maps between day and night at 12 or 36 hours 
intervals . The HCMR experiment was originally designed to produce thermal 
inertia data for soil and geology applications but the very good performances 
of HCMR are suitable also for oceanographic studies. Data were received from 
NASA (National Administration for Space Research) through an investigation 
concerned with sea surface temperatures of the coastal zones of France. 

Available HCMR data are photographic or digital products covering 
a 700 x 700 km 2 square scene. *£he following informations are displayed : 

(1) surface diffuse albedo or reflectance in the visible channel (day only), 

(2) surface temperature from the infrared channel, £3) surface temperature 
difference between day and night, (4) thermal inertia, which was not used 

in the present study. About 1000 scenes covering the coastal zones of France 
were received for the period May 1978 - May 1979. Examples of the photographic 
products are given for two areas in the Western Mediterranean Sea (Fig. 1) in 
the North Sea (Fig. 2)) where large diurnal variations of the SST were ob- 
served . 

II- 2 - Diurnal heating and glitter (sun glint) patterns 

A large number of the received date from May to July 1978, over the 
Mediterranean Sea exhibited very interesting and similar features in both 

A 

the visible and the infrared channels, as shown in Fig. . between Corsica 
Island and the south coast of France, and also close to the east coasts of 


of Corsica and Sardinia Islands. 


Warmer araas in the thermal chanel are associated with changes of 
brightness in the visible. 


The observed changes of brightness in the visible are Identified 
as glitter or sunglint patterns - i.e. specular reflexion of direct solar 
radiation by the wavy sea surface. During the concerned period around the 
summer solstice, the observation angle of the HCMR imagery was allowed to 

* 

be very close to the specular reflection of direct solar radiation, in the 
western part of the scenes, which is favorable for 'observations of glitter 
patterns. Most of the time, the glitter increases from rough to calm seas, 
when the wind decreases and the sea surface becomes more specular, and exhi- 
bits a maximum brightness when the observation angle is close^he specular 
reflexion of solar radiation : a homogeneous bright area is thus noted in the 
south-west part of Fig. 2a. But for very calm seas, the surface reflexion 
becomes nearly specular, and a decrease of the brightness may also be observed 
because it is very unlikely that the observation angle is strictly towards 
the specular reflexion. Such a darkening is observed in the northwest part 
of Fig t J2&, where the two processes are present with both bright and dark 
areas corresponding respectively to weak and nul wind speeds. The fact that 
smoothing of the surface could produce either an increase or a decrease of 

&ztU. 


the glitter brightness was previously mentionned by LA VIOLETTE asMMHh 
(1980) . A physical and detailed description is given in Appendix, to support 
a further quantitative analysis of the data. The dark patterns in a mean bright 
glitter can thus be clearly interpreted as nul wind and calm sea areas, which 
obviously are favourable to a larger diurnal heating of the upper layer of 
the ocean because the heat transfer to deeper ocean layers is limited by a 
low turbulent mixing and thermal diffusivity. 
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II- 3 - Meteorological observations 

Evidence of a large diurnal heating corresponding to low wind speed 
conditions is also given by correlative meteorological observations, Surface 

observations are presented in Pig. l-.ii for the case in the Mediterranean Sea, 

O 

and in Fig. 2-d for an other case found in the North Sea where, due to higher 

A* 

latitudes, glitter is almost always unobservable. On Fig. 2-^ a large warm 

l 

spot was detected by HCMR in the middle of the North Sea wich is coincident 
with the center of high anticyclonic situation pressure whertnul wind speed 
is reported. Warmer areas observed in the Mediterranean Sea on Fig. 1-b are 
also coincident with low or nul wind speeds, but the observed wind field is 
much more complicated because most of the reporting coastal weather stations 
are affected by some breeze effect*; wich surimpose to an anticyclonic circu- 
lation. Cloudfree satellite SS? observations are frequently acquired during 
similar anticyclonic situations with moderate wind speeds, It must be out- 
tined that satellite estimations of SST may thus be systematically affected 
by diurnal heating, and a tentative statement of this is discussed in 

section 4* 

11-4 - Day-Night observations 


Heat loss during the night very rapid^ly destroys most of the diurnal 
heating, at least in the upper layer, which wa6 produced during day time. 
Evidence of a diurnal heating may thus be found from a comparative analysis 
of two successive day and night observations at 12 hours intervals. For the 
two cases given in Fig. 1-c and 2 -b, nightime observations show a much more 
constant SST field and the noticeable daytime warmer features disapenr. 


Figure 1-d gives the result of the computed day-night temperature 
differences after the proper calibration algorithms have been applied by NASA 

. These differences present the advantage to be independent 

of the mean mesoscale SST field and allow to enhance the diurnal heating, which 
again closely correlate with glitter patterns in the visible channel. Day-night 
temperature differences are used in the followings for a more quantitative ana- 
lysis of diurnal heating. 

» 

III - DEPENDENCE OF DIURNAL HEATING ON SEA STATE AND WIND SPEED 

I 

The observed diurnal heatings were further quantitatively analysed 
to derive its relationship With the sea state and the wind speed. Day-night 
temperature difference were correlated to the reflectance of the 0.5 - 1.1 uo 
channel. This reflectance, mostly due to sun glitter, is related to the sur- 
face slope variance and to a mean wind speed using the statistical model from 
COX and MUNK (1955) . 

III-l - Diurnal heating and glitter reflectance 

Day-night temperature differences (Fig. 1-d) - i.e SST diurnal vari- 
ations - show patterns Similar to the glitter patterns (Fig. 1-Jrf) , on June 3, 
1978. Fig. 3 gives the result of the correlation obtained when the diurnal he- » 

ating, AT, is plotted as function of the glitter reflectance, f) , in a small 

« 4 

study area, east of Sardina. Most evidently a close correlation exists and AT 
rapidly decreases when p increases. To further interpret that fact, p. has 
to be related to the wind speed, or more exactly to the statistics of surface 


slopes . 


Using the statistical distribution of surface slopes from COX and 
MUNK (1955) , a model was developped to relate the glitter reflectance to the 
wind speed. This model is detailed in Appendix. Results indicate that p could 

either increase or decrease with wind speed ; p presents a maximum value 

J s 

for a given wind speed value which both of them depend on solar and observation 
angles through © n (tg © n is the surface slope allowing specular reflection to- 
ward the sensor) . Fig. 4 give the relationship between pg and the wind speed, 

U, for © n ■ 8°, 10°, and 12°, which correspond to the area previously studied 
for AT » f (p^) . In this case p g increases rapidly at the lower wind 

speeds and then is rather constant for U > 3 m.s” 1 so that U can be estimated 
with a good accuracy from p^, only for U < 3 m.s”*. The study has thus to be 

limited to this wind speed range. It should also be noted that p is physically 

9 

linked to the surface slope variances, and only statistically to the Wind speed. 
Local anomalies may thus occur, in particular when the fetch of the wind over 
the sea is variable . Keeping in mind these cautions, we may now transform 
AT(pg) in AT(U) which is given in Fig. 5. 

I I I- 2 - Diurnal heating and the wind speed 

The first point to be noted on Fig. 5 which gives the diurnal heating 
as a function of the wind speed, is that AT rapidly decreases from several c C 
to 1° C when U increases up to 2 m.s *. The scatter of observations on Fig. 5 
is remarkably less than on Fig. 3 for AT(p^) , because the variations of 

p g with changes of observation angles within the study area have been t iminated. 
A fit of AT (U) on Fig. 5 would give : 

AT - 0.4 U + 1.1 


(in °C for U in m.s *) 


( 2 ) 




J. - B - 


Some uncertainties related r,o the model p (U) have been previously outlined. 

V 

Additional errors may be due to atmospheric effects on the measured radiances. 

An aerosol atmospheric reflectance§4bout 0,02 was estimated from the minimum 
reflectance within the scene (p g « O) and substraqted in the 0.5 - 1.1 pm 
channel. Day-night temperature differences have not been corrected for atmos- 
pheric emission in the infrared. This approximation would be valid only if the 
atmosphere remains the same between the two satellite overpasses, but a bias 
due to a change of atmospheric parameters - i.e temperature and water vapor 
concentration - could have occur which would possibly explain the 1.1° C cons- 

* 

tant found in (2) . Last, the observed AT are certainly underestimated by a 
factor t» the atmospheric transmittance in the 10.5 - 12.5 pm, which is typical- 
ly x* 0.7 for a midlatitude summer atmosphere. 

The results may be compared to the predicted values from HASSE (1971). 

_2 

Using a mean solar irradiance at Je.3 level Q ■ 900 W.m in (1) , AT is found 

n 

^ — 1 •* i 

to vary like U (U in m.s ) which fits the measured values in the wind 
speed range 1-3. m.s"’ 1 , but overestimates AT for U < lm.s 1 . As pointed out by 
HASSE, the results of the model given in (1) can not be applied to the lower 
wind speed range because the model used by HASSE refers to a steady state as- 

<oi 

sumption which is then not respected at scale* of a few hours. 



Fig. 5 and other HCMM scenes with large diurnal heatings indicate 
that diurnal heating do not exceed about 5° C, and that a limit value should 

exist at low wind speed. This value may be obtained by solving the heat transfer 
equation : 

(ku> = pc (3) 

dz dt dz dt 


for k(z) ■ k the thermal molecular diffusivity - i.e no turbulent diffusivity 
at U * o. Eq, (3) was solved using the following conditions s 

F ( 2 ft) - F (o , t) g (z) - F Q (4) 

where F(o,t) is the solar irradiance at sea level, F q the heat loss by the 
surface, and 


g(z) ■= l a i axp(-k.z) 


(5) 


where a^, k^ are given in Table 1 and were obtained from a fit of g(z) accor- 
ding to the work by PRUVOST (1976) . g(z) is taken as independent of time 
in (4) which is a rather good approximation since the underwater penetration 
of the direct solar radiation is close to the nadir even at low solar elevation 

angles. An homogeneous layer defined by F(o,t) g(z o ) ■ F q was set just below 

dT -2 

the surface for which (•?—) * o (z„ is a few centimeters for F * 100 W.m , 

az o o 

-2 

F (o, t) » 1000 W.m ). Under these conditions, AT was found to vary nearly 
with the net heat budget of the surface : 


AT a C 
max 


(F (o, t) - F o ) dt 


( 6 ) 


—6 —1 

with C ■ 0.65.10 K.j m . For the HCMM observations or 


✓ 




f: 


! 1978, 


-2 


(F(o,t)-F o ). dt was estimated to about 600 W.m during 4 hours (in fact 

-2 


a maximum value of 900 W.m at noon at satellite overpass) and 


AT 


max 


■ 5.6 °C 


(7) 


The Hasse's formula (1) may ba simply accomodated to account for the 
limit found in (6) by writing : 

AT - 3.5 lo" 3 Q / (U it,) + U) (8) 

where U Q (.tj wi dependtof the given hour during the day. In our caie, U Q should 
-1 

be about O.J- m.s and when plotted in Fig. 5, Eq. (8) fits pretty well the 
observations . , 

* 

III-4 - Frequency of diurnal heating 

From May 13 to August 28, 1978, 60 HCMM scenes taken over the Westernn 
Mediterranean Sea were examined, of which about 34 scenes exhibited large diur- 
nal (-typically more than 1° c) heating of particular areas of 10 to 100 km width. 
Relative frequency of the event is rather large, and is enhanced in some areas 

i / 

affected by a bree.ze effect where the wind systematically becomes nul at some 
distance of the coast. Table 2 give relative frequencies of low wind speeds 
(U < 3 m.s ■ ) at some stations along the Coast of France during the summer 
months (from DARCHEN (1974)). Frequency of nul wind allowing a diurnal heating 

of more than 1° C are between 10 to 30 % .Frequency of low wind speed 

_1 ' 

(1 < 0 < 3 m.s ) is from 20 to 50 %, allowing a diurnal heating of about 1° C. 

From these frequencies, and N^, a mean diurnal heating AT was calculated as . 
AT » 2.5 N 1 + N 2 

and is given also in Table 2. The mean diurnal heating ranged from 0.5 to 1.5° C 
along the south coast of France with a maximum on the French Riviera (Cap Ferrat) 


IV - CONCLUSION 


The present investigation, using SST satellite observations from the 
HCMM experiment has shown a high frequency of large diurnal heatings (more 
than 1° C) of the sea surface during summer months in such areas like the 
Mediterranean Sea where low wind speed are very frequent. This shows 

that satellite observations at noon and during the oftemoon should be rejec- 
ted, or at least checked to eliminate those corresponding to low wind speed 
(U < 3 m.s S . If not, a systematic bias could be introduced in the SST ana- 

lysis of s ^as areas, particularly the marine coastal areas affected by a 

1 } 

sea-pand breeze effect. 

Using simultaneous observations of the glitter reflectance, the diur- 
nal heating was correlated to the wind speed. Diurnal heatings of about 1° C 
were found for U = 2 m.s* 1 , which fits the formulation given by BASSE (1971) . 

A maximum diurnal heating of 5° C is found for nul wind conditions, which is in 
agreement to the value calculated from the radiative and heat 

transfer, assuming the thermal diffusivity is only molecular. 


APPENDIX 


, Glitter refers to direct solar radiation reflected 

by the sea surface. This reflection is specular for a planar sur- 
face. When there is wind, the surface is agitated and consists 

of elements which are statistically distributed around the horizon- 
tal plane. This produces a more or less bright spot of variable 
dimensions which is commonly called glitter. 


The radiance L reflected by the agitated sea surface 
can be expressed (COX and MUNK, 1956) 


L 9- 


E g R(u) 
4 w v^n 


and the equivalent reflectance 
ILL 

_ H R(to) 


y , E s 


4 u u u 
M s M v H n 
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Pg will be expressed as 


(A-l) 


(A-2) 


where E is the direct solar radiation at sea level, 
s 

R(w) is the reflection coefficient of water at a given indicence u>, 


p is the probability of encountering a properly oriented surface 
element , 

y - cos0 / y - cos0 , u “ cosG , respectively define the zeni- 

V vs s n ft 

thal angles of the observation direction, the direction of incidence, 


and their bisector, 


v> is the angle between the planes of incidence and observation : 
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CO 3 ^ 
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Prom a study of aerial photographs of glitter patterns, COX and 
MONK (1955) developped d in a Gram Charlier series which in a first 
approximation is reduced to a gauss i an distribution , with re- 


volution symmetry : 

(W 9„) 2 

— Z“ 

Ha" a 


p ■ — iy exp - 


(A-5) 


with a 2 - 0,003 + 5,12.10~V -1 ± 0,004 

m.a 


(A-6) 


for l<y<i4 m»s 
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Figure 6 gives an example of the glitter spot Og 
thus computed as a function of solar zenithal angle for different 


values of W, and for a nadir viewing (© v ■ 0) . In accordance with the 

reciprocity principle, by permutation > Fig. 6 also gives p^ 

as a function of observation angle for a sun at the zenith (0 - 0) . 

For a given angle p presents a maximum, o , at a certain value of 
■ g ■ gm 

a which is related to wind speed, a and o are given by : 
m in gin 
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ha dashed curve in Pig. 6 envelops the proceeding curvec and repre- 


sents the maximum glitter p^ as defined by (A-8) . 
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Table 2 - Relative frequencies of low wind speeds : 

Nj ! nul ; Nj : Beaufort forces 1 and 2 (1 < 0 < 3 m.s 

during June, July and August in the french mediterranean 
coastal area from DAJRCHEN «An estimate of the mean diurnal 
heating St is given in column (3) . 


Station 

N t % 

N 2 % 

St 9 c 

Cap Bear 

16.0 

26.9 

0.67 

sate 

9.5 

42.3 

0.66 

Panagues 

21.3 

26.8 

0.80 

Cap Camarat 

10.3 

46 . 6 

0.74 

Cap Perrat 

35.1 

50.4 

1.38 

Cap Corse 

18.4 

35.5 

0.82 

Pertusato 

6.4 

21.0 

0.37 

; 42° N-6E 

; 7.6 

1 

i. 

/ 

1 

0.5 ? 



FIGURES CAPTIONS 


Figure 1 - Diurnal heating in the Western Mediterranean Sea j 

(a) 1 HCMM scene A-A0038 - 12440 on June 3, 1978 at 12.40 TV, image 

w 

center is at 40.54 N, Oil .04 E, Visible channel s darker tones 
are is lower reflectances. Hots the bright patterns East and 
West of Corsica and Sardine. 

(b) - Same as (a) but infrared channel : darker tones are colder 

*W 

sea surface temperatures. Note warmer waters East of Corsica 
and Sardine. 

<c) - Day-night temperature differences between HCMM scenes obtained 
on June 3, 1978 at 1.50 TU (night) and 12.40 TU (day) . Darker 
tones are smaller diurnal heatings. 

(d) - Meteorological situation, on June 3, 1978 at 12.00 TU. 

Figure 2 - Diurnal heating in the North Sea : 

(a) - Day HCMM scene A-A0034 - 13120, on May 30, 1978 at 13.10 TU. 

Image center is at 54.27 N, 00.01E. Infrared channel : darker 
tones are colder waters. Note the warm (bright) spot between 
Scotland and the top right of the image where a thermal front 
is shown close to Norway. 

(b) - Night HCMM scene A-A0035 - 02280, on May 31, 1978 at 2.30 TU. 

Image center is at 56.13 - 03.00E. Infrared channel : darker 
tones are colder waters. The warm spot disappeared during the 
night . 

(c) - Meteorological situation on May 30, 1978, 


Figure 3 - Day-night temperature difference vs glitter reflectance on June 3 , 
1978 , for a study area East of Sardina, 

Figure 4 - Retrieved wind speed vs glitter reflectance for the study area. 

Figure S - Day-night temperature difference vs retrieved wind speed for the 
study area. Dusked line is from MASSE (1971) . Full line is (8) i 
the KASSE's formula after modification to account for a low wind 
speed limit of AT. 

Figure 6 *• Glitter reflectance vs zenithal viewing angle, for a sun at zenith, 
and several wind speeds from 0 to 15 m.s"’*'. Maximum glitter reflec- 
tance is given by a dashed line. 
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APPENDIX D 


The QolZouUng tuning give ike date, Identification and location of 
eenten of Image of HCMM a tenet received from NASA by the Pntnclpal Invent- 
gaton. The loti column "ETAT" give ike ttoau of ike connet ponding digital 
data : 


- R .• received 

- IR .• received but not readable 

- C .* nequetted but not received. 
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7DEC75 
-7.de:? 3 
7 d E C ** 5 


??4-l3440-2 

22«*-i3*$o-2 

225-1 2230-2. 
_ JLr r '-* 2150^1 
223rl223 0-V 


'230 -2. 


225-12250-2 

225=J2260-JL 

225-12260-7 

223-1 2270- 4 ' 

225-12270-2 

225-52280-1 

225-1 2280-2 


a3.20n 

i?«a_Sn 

. 37; 5 5n 
*3.1 1*..-. 

* l*ain 

4i.ii . 

*6_»05n .. 

*»5.0'5n 

4^l4*_ 


2 *• . f) 7 1, 


'6DFC-3 

1 fc D F r - s 

' 7DFT-? 

' 


• 13370-1 


-*• 7- 


2* ?*»- 
r * t 


11 . 02 v 
13 . 1 Ow 

11, 10W 

-1&.12E 

.30. 32* 

p . : c r 

Bv442 

5.44*. 

7.C2E 

7.02E 
6.37E 
-6. 37E 
. _4.34E 
4.14 E 

6.. 06F 

.fcvOoE 

3 . 4 6 E_ 

_ ^46t 


‘.7. 5 IN 

*7 .58*; 

>4. DON 
P4.03N 


» ?• . ? ’ v 
* . ? 7 N 


: / 

r 


\in 


r . 4 ' E 




$ r E i> E 


z A 1 1 


iDEKTJf JC4TID5 


LOCaTICN 


BCE ETfcT DST Pk 


'.7{J Er**S 

235-12132-2 

. *1 5 

-9.49E 



1?0EC78 

235-12150-1 

54.-2N 

.7-17E 


^ • i. r ■ 

17PEC73 

235-12150-2 

54.42N 

7. 17E v . . 



i £0 E C78_ 

236-12260-3- 

— 9„JCE 



1 6D E C 78 

236-1 2 2? 0-2 

36„24K 

_ 9.10E 

jj— 


.1EC E075 

_236-l23l 0-1 

*5.325 

$+ 1 9 1 



1 SPECKS 

236-123*0-2- 

4* 8.32V 

5. 1 9 E 



?80EC*»3 

236-12333-1 

&4*34N 

_ .. _ 



‘•SOEC-'S 

236-1 273''-: 

>4.T4N 

7.47E 



i9X>£073 

23’-124aO-1 

.36. 07N 

6*41.1 


, . — . . 

* 9p£ C7S 

’7?-l2*.6P-7 

. o?n 

4 . 4 1 £ 


• 

i9*e:-3 

22"-1 2h£0-1 

-:.i2'. 

7 .55 £ . 



'9DE078 

?y*-l24S3-2 

- 2 . * : n 

?.55£ 



' ©oE:?i 

.237-12490-1 

-s. 1 6.; 

-521 

e 

w 


^Coe;*»3 

2’*’-1 2490-? 

48,14m 

,52£ 

c 


i ODE :?s 

237-12510-1 

)4,14N 

1.26* . 

.X 


''99E0‘’S 

2^7-12510-2 

7- . 1 C>5 

. 1 * 3 6 W . .. 

c 

. • 

2f'DEr75 

. 238- *>3£0-2 

37_23I 

135. 0-E. . _ 



21 DE075 

27°- 2290-3 


7.5*W 



? - d £ : *» 3 

its. ? 2 ^ 7 

43.495 

e.004 



2ide:?5 

23 5 -l 3250- 1 

-5.055 

7.04*' 



2ixe:?s 

230-12250-2 

•o. oSn 

_ ? « C 4 * . ... ... . 



2 2 0 1 0 ' 8 


*9. 77s 

17. 16* 



?2?e:v» 

24.0-13—0-? 

-*9 , i7:. 

V* . 1 ©* 



?2?s :?s 

?41-12?70-1 

-2 .74\ 

- 5 .5 5 E 



22D E C?i 

. 241-12230-2 

42.34V 

®. 55E. 



? 4. 0 E 9 •» * 

242-124*' o-i 

39,< *f; 

5 . ??E 



752s7*M 

Z- 7 - 20-0-7 

- 4 , 1 9 1\ 

1.4?* 301 

R 

20 

:8u:'-3 


* 9 , % 7 ri 

. 77 E 



7soe:7s 

2-3-12530-2 

39.5 7’; 

. 77£ 



250EC-3 

247-130*' ''-l 


i . 1 ? u 3 01 

R 

k 14 

25DEC73 

247-13000-2 

46,31 • 

*.18* 301 

R 

k 17 

250EC73 

247-* 3120-1 

52,36* 

3.36*' 



75or:73 

2- "'-l 3020-2 


7 . 3 e *' 



?6se:7S 

?44- 2210-7 

-5 . ?4\ 

5.5 3*' 3 01 

R 

H 2 3 

?"of :**s 

. 245- 2400-’ 

-4.00*. 

1 9 . 5 5 w 



2Et>E C78 

24*- 1210-7 

<•0.525 

7.52E 

C 


2iD£ C*»3 

. 246-12160-1 

- 3 , 3 1 :> 

10.11E . 

.c 


zsocCrs 

246-12160-2 

-3. ’-If; 

10.HE 

r 

w 


:®de:* , 3 

24'- 1-00-7 

37.13M 

2.14E 



790EC-3 

247-", 2360-7 

-8.4 35' 

■» . 51E 



2 C ?E 0 " 3 

2-'- l27’0-' t 

9 - .4 - 5 

1 ,19E 



3 f! 0 E 0 i* 8 

2- p -1 25O0-* 

35.555 

3 . 2 OE 



3?0S:78 

24.7-1 '^20-1 

42.015 

1 .36E 


f 

" r* D E 0 ■* ? 

? 4 * - 5? 0-2 

- 7 . r\ 1 \ 

* . 36 E 



3* 5E — S 

^ Z ^ 

4. «• 

45,53V 

4.07* 



T 0 ; i [ " " • 

24 0- 

- e . 4 V. 

4 . ?7» 



v* v- ir-2 

2*3-1 3i 7 >> 7 

j* V: 

6.2-* 



^ i * ?f) 

^ J * • 7 

n?^- 7 

-1 .73’. 

5 . T 7E 



7,ji .'-9 

* 7 1 * 7 

4. ^ . • • »*l , * 

35 . 775 

' .III 




r 


CATE 

identification 

LOCAl 

3,-AN79 
-- _lJ M?9_ 
. ^ .3J*N'79_ 
3 J AV7l_ 

252-12260-1 39j56N . 

23^12210-2 _ i 9 .56'v - 
252-12280-1 *6 .t:n 

2 5 2- 1 22 8 Q --2 .4 6- Ml 

Zj A 

_3j am**9 

252-12290-1 
252-1 £290-0 

52.-0?* - 
. 3’^.P?N._ 

t jak-»? 

Xj an ?5 

55- 2260-7 

.-62-. 50* — . 


- - . r 7 •'• 

A JAK79 

255-12210-? 

LU ^07k. 

7 j A u 7 9 

._256-_1060-3_ 

i2J*l 

"7 JAN79 

..25«- 2410-7 

S ? . 4 4 N 

_7 j A N 7 5 

_2 5.6- 13395-1 . 

^3.52n__ 

7.JAM79 

25 6-J 3390-2 

- 2 . 0 2 ■ • 

.. .. 10JAN79. 

_ 25 5- ----- 

32.19N 

1 C J AK79 

255- 2 4 -o-3 

il . 4 ON 

51 J AN*7? 

.. 260- 2150-2 

S3 .z 4 ^ 

1 2 J A N! 7 9 

261- .2220-2 

3 m . 4 2 N _ 

_J»3j AK?£, 

_ 26 2-. 11* 0-2 . 

. *• 1 ,2iJi . _ 

13 JAN79 

: , 2- 44 c^7 

3 5 . 1 9 n 

1 3 J A N' 7 9 

262- 2520-2 

: 0 . T * ’■ 

14JAN7 9 

263- 1220-3 

3?.r5N 

Uj AN 7 ? 

_2t'»- 1310-2 

39.55V. 

1SJAN79 

264- 1500-2 

53.47m 

15J AK?9 

264- 1510-3 

- 7 ..AAN 

^_i.S J AN 7 9 

264- 1 570-2 

* 1 — "7 9 N 

r 1 6 J A M - 9 

265- 22.90-3- 


1 6 J A N 7 ® 

265- 2110-3 

42.7CN 

4 7 J A •: •» 5 

2t6- 2260-2 

3 2 . 4 ON 

1 7 j A \ 7 9 

266- 22*0-3 

*• 6 . 7 7 N 

4 ?J A 

266-1 2220-1. 

40.42V 

17 j AN? 9 

266-13230-2 

40.42N 

<7 J AN79 

266-13250-1 

4 6 . 4 6 

17 J AV76 

266-13250-2 

46 .46N 

4 7j A’.'-9 

266-1 3 260- 4 

3 2 . 4 9n 

17 J AN7? 

266-13260-2 

92.49N 

4 £ j « N 7 c 

26?-l34?0-< 

-2.55m 

1SJAN'? P / 

26?-1 34?0-2 

-2.55m 

1 9 J A \' ? 9 

_2* p - 1290-2 

35 . 74n _ 

21 JAW79 

2?0-...-20i Or? 

6l^i?N 

2 1 J A N 7 9 

27?- 2 01 Or 1 

6 4 , i £ N 

2 4 J A N *» 5 

2 7 * - j * t n «. t 

4 5 . 1 4 N 

2 1 J A \' 7 s 

27?- 2340-3 

39.05M 

21 J A\'75 

27P-13900-1 

46. ION 

2 1 J A N 7 ? 

27*-1 30C?*7 

46 . 1 ON 

2 4 J A N 7 5 


S2.12M 

2 1 j a *; 7 * 

2 # - 1 3 * * 2 

92.1 ? N 

?2 j t : m ~ 

2^2-12*'? j-1 

51 . 1 3n 

* 7 t 


4 * « ^ 

+ % — 

* *7 • • 

7 * J * « y 

j - *--** 7 * 4 "-* 

-4 .7 9*. 

2 2 J * N 7 * 

C < «. »•*•)•> — 

-4 . 3 VN 


scene 


SDE ETAT DST P* 


f.25t ... 




6^J2.9E_ 


.4-3TS 

ti t lJ 

6.27A 


6 .*7* 

i .-w 

.1 0-.3&H.. 

1 C.'SBW 

2^3. 1i_ 
- . 55 * 
. 1 - 4.9 4 __ 
_5.4?v . 


3.01_...l r 49 


.7^.041 
12.2 5 w 
8 ,1?E 
_3.S5E 
X.25T 
i . 5 7 r 
.CCI 
_ 2^ 1-5 W_ 
va.2Dv’ 
L. 57.' 
7 . 20 * 
5.56* 
5,56* 

7.5 5W 

?-» 5 5 W 

1P. - «N_ 

If , 18w 
1 4 . 08w 
1* . ORw 
.4.32F 
_- .-36E 
. 36E 
1 .4 OK 
7.35* 

1 .364 
4 , 3 6 * 
t . 5 6W 
7,5ft* 
V . *5E 

* 4 * e r 

» » ^ * 

* * .-?* 

4 * . 47; 


302 
2 0 2 


K 

M 


2 

5 


_C 

* 


L 


<jy 


CATE identification 


LOCATION 


Sf ENE 


bCE CTAT 1>ST P* 


23JAN79 

272-13370-1 

92.43N 

u.eitt 


23j» v?? 

272-1 3370-2 

?2.43:. 

-13-01W . .. _ 


2m J AN 79 

273- 1210-3 

* 1 . 20 * 

7-01E 

c 

!4JiH7t 

273- 1220*3 

. 35. 1 3n 

5 - ItE _ 

T .Z ""C * 

25JAN79 

274-12340-2 

-40. 3** 

5 . 33E 


25jAN'*9 

ZiJAN^C 

, 274-12360-1 
_27a-123cO-2 

. 46.44N 
*6.4-K 

3-34C 


2.34E 


25 J > N79 

2 *--i 2 ??o-i 

> 2 - 47m 

1 . 1 1 i 


25j AN7« 

274-12370-2 

52.4?a 

T. 11 I 


2tJ AN75 

. . 275- 1 590-3 

38.38* 

2-574. 


2 6 J A N 7 ® 

275-1 25?0-1 

5 1 . 4 OK 

1 .UE 


U ji'.s 

275-12570-? 

J5.40K 

t-UE . 


26JAW79 

275-12550-1 


7.564 


_ 26jA\7® 

275-12550-2 

51 • 4 9N 

2.364 


2?JA\'79 

2 ? 6 - 2 U 0 -? 

>C.4*V 

3. 294 


_ 2£jA,'7® 

27*- 2330-2 

-7. 

.. _9-l6W 


28 J AN'70 

2’ 7-13300-1 

•44 . ?9n 

9.264 


23JAN79 

2*7-12300-2 

4-.?9N 

9-264 


28 j a 

277-1 33?0-1 

92.33k 

11 .394 


28j A ."9 

777-13320-? 

5*. 33* 

* • ,394 


7-5 j * .-9 

273- 1170-3 

36.1 2N 

6.511 


3 2 J A *?9 

279- 1320-3 

-3. 00n 

. -4.16E 

. c 

3 ? j A S' 7 9 

270- 1355-7 

36.540 

? . 2 8 L 

* 

% 

3 * j A v * 9 

279-12230-1 

* 9 . *•* 

6 . 5 1 E 


3D J A NT 9 

270-12230-3 

4%?S* 

fc.5iS 


31 J AN' T 9 

280- 1490-2 

7 1-77* 

2-491 _ . 

c 

? * J A 7 9 

28*- 1510-3 

*5 . S3* 

.31 E 

* 

V 

3 « j a y 9 

23 0- 1 5 *» 0 - *» 


1.23'/ 

c 

1 F£?’i 

?«*- 20S0-3 

51 . ,' 6 \ 

1.5*4 


1 FES?® 

201 - 2*9 0-’ 

-5.?3V 

4.27- 


2 F e 9 7 ® 

?« 2 - i i - *» 

-9 . 1 5N 

7.194 



283- 1*90-3 

*•? , 4?N 

1 0 . CEE 


3?E?**V 

2 o’- 1110-3 

36. 4 IN' 

« . 2 0 E 

c 


254- 1 2 5 0 - 1 

> 0 . ? 6 N 

* . 1 6 fc 


-foo*? 

25-- 1270-? 

-4.22N 

6. CAE 

c 

4 F £?’“• 

25-- i’oC-3 

35 . 1 6N 

4. HE 

r 

w 

* r S ■»*»9 

2 34-1 2220-1 

39. 5 IW 

«.A«F 

c 

,tEi?c 

234-12720-7 

3® ^ *1 o 

8.41E 

m 

W 

5*0*79 

?p«- l-’O-? 

7 2. < ON 

3 . 5 0 E 


5 r r « 

2? 5- U50-? 

-4,46V 

1 .37C 


5 ► E 0 7® 

755 . 1475 — 3 * 

3F.41N 

, 1 6 F 


3FES-»® 

2*5-12450-1 

35. niN 

4 . 2 4 E 

C 

5FE°7® 

275-1 2400-2 

33 . *1 N 

4. 248 

r 

% 

5 F E 9 79 

235-12430-1 

7 * . 1 1 V 

. 22E . 


5 Ft?? 9 

7ce.i 2475-2 

?0.«*V 

.221 


f rr**' 

r?A- 7*1*. - " 

5 1 . 4 * V 

. 2SC 

* 

ff£ V ’' 

:#*- ? * ? 0-7 

1 

* .5®. 

m 

^ • 7 5 •» ,' 

:?7-i?57:-i 


. ? r t 


a r?^; 

• • w • £ . 4 

4 / * *w ,» 

. 2or 


7rf a-*; 


3*. 77k 

5.134 






date IDENTIFICATION lCCaTicn 


257-13160-2 
_ 237.-131*0-1- 
IITrUUO-2 

_ 290- 1170-3 
2«4-12**0-7 

.. lUsl IJtfisJL 

?‘?5rl2260-- 
2 o?-« :?t0-2 
96-1 


39,?7n 5.13w 

* 5 - A 3n 7 ,06* 

45,c3n _-^,Of.w 
S-JLT. 

A2,f>4* 11.7AT 

A?_.p4N 1.L.14E 

35.A6K 7 . *• ' 1 

3S..X6N L.AJ&I 

* 


^7FEe?9_ _ 

297-1302Q-1. 

.3**04k 

*7f E*79 

297-1322^-2 

36. oak 

l7rjo*»9 

297-130*0-1 . 

a 4 .» 1 0.N 

«?PEB79 

2«7-i 3%0-2 

44.10* 

2:FE3?_9 . 

213-1 2103-1 

-38..4 9N 

’0FEB79 

203-1 2*90-? 

36.49*. 

2Z F £ 8 ?9 

320-1 £210-1 - 

- - . 5 5 H 

2PFFB79 

300-1 2?^ 0-_2 . . 

.‘•4.55N 

?rr.£979 

300-1 2220-1 

70.5°4 

20FE979 

?1 c E ? 7 v 

220-1 2 270-2. 

331 - 1 27 f 0-1 

70.59K 
a 1 . 7 5 V 

?l FE379 

221-12^90-2 

*. 1 . r 5 K 

71**379 

3«i-1 ZAr^r* 1 

a? # a0N 

?1 CE^O 

?2*-1 2a20-2 

a7,40K 

22FE«*»V 

322-1 25* 0-« 

3 ? . 5 3 N 

22f£P75 . 

302-12590-2 

36-535 

22FE9J9 . 

.202-1 2590- 1 _ 

7.L..04K 

72f E‘'* , 9 
PlFt"*?? 

? 0 ?-l 25 e?»“? 
3* , »* 31 

5 ^ , 04 K 

-A . • ‘ f 

2 3 f E ? ' 5 

t» r f >.5 

aa . 1 - \ 

23 F £ a ^ 5 

3*!2-l3 1 6 0-i 

30-lfcK 

23 f E 3 79 

302-13160-2 

70.-6K 

2iFE97? 

37A-1 1550-1 

32. A 25 

24FE97V 

?34-115 f 0-7 

37.42k 

2 •♦ r c 3 ? V 

2aFE«>*»9 

306-1 37j?-i 
204-13320-2 

a* , *9\ 
44,795 

2 5 F EE79 

305- 12130- * 

A 2 , A 6K 

25FEF79 

225-1 2120-2 

A2.46K 

2 5 f F ? 7 f 

325-1 21 5 0- 1 

-5 ..5 2*. 

25FEB77 

225-1 2150-’ . 

A.*. *7* 

25 F c 9 7 9 

205-1 2160-'’ 

54.5AN 

25FEB7V 

2^5-12*60-? 

3A.54K 

?6CE 6?7 

? 06-1 2? r 0-1 

3 0 . 1 A K 

7 6 2 E 3 7 9 

306-1 2 ?, 00 -? 

3« . 1 4 It! 

< C •) 

? £ F c " * 

“ 6 r E 3 ^ ' 

? c r ? a - ; 

5 pr* 1 ( 

366-1 2370-7 

* * a- 1 77 ? 0 - 4 

4 3 • 7 

4 5.ro:.' 

3 * . 7 A ‘ 






: A T £ 


KEMIFICaTIOK 


SCENE 


EDE ETAT CST Pk 


tCCtTjCN 


27(E*?9 

•<77-12480-2 

32 . 01 * 

2. OoE 


27FE?*»C 

T-07-12%90-1 

-4.*7K 

- .,ni 

c 

27tES79 

377-12400-2 

44„fl?« 

.161 

.. c 

27?E»7? 

377-i25lAf3_ 

)LJ1I 

1-56 * 

... . S-. ... ... 

27FEP70 

3j7-l25*D-2-_ 

>0.11* 

-.56* 

z 

1**A®?9 

_ 30O-13250-* 

_* 5..11 N 

LMI 

mt • «• • a. , ■ , _ # ^ 

-1v AR?9 

309-1 325C-2 

■*3. Pi* 

9. 03* 


1 4* AR*»9 

_ 312-12^70-1. 

3?.4lK 

6*07 1 

#» 

... 

1 410*^9 

322-" 1 2270-2 

3814?* 

4 . 07E 

c 

« 4 *» A 7^0 

322-12290-2 . 

*4. *7k 

- . 1 3 E _ 

c 

* $m A9*9 

223-1 2*e0-”i 

39.55k 

1.0SE 


1 5**17? 

323-1 2460--2 

J9.55N 

1.081 


*6**179 

3?a-i3040-i 

39. 7 4 K 

3. 1 1 W 


16* A0*»P 

2?4-i 3040-2 

39. 7*N 

3.11* . . _ 



7?4-l3»)50-l 

•*5. r9x 

5.05** 

c 

1 6* A $ ?9 

324-13050-2 

4,3. r9K. 

. .5.05*’ 

c 

1 7«* aR?9 

325-13270-1 

* 1 .rlK 

P..20* 


* 7*A1*9 

J23-i3r7o-? 

* 1 . r> 2 * 

1 .2Cr 


1 7mAC79 

325-13240-* 

*7 , 0? K 

1C.20W 


1 7v AR’9 

325-1 32-0-2 

47 . 

10. 20V 


1 $ u A e ? ? 

326-120*0-1 

43 . 7 H 

* r . 35E 

r 

V 

1 fiv a»79 

326-120..O-2 

43«2tl 

V.35E . 

_.c 

* 1**179 

326-1 2060-* 

55.74U 

5 . 5 0£ 


1 £*AR?9 

326-l2r.f.o-? 

. *-?• 

5.5CE 


1 Q* A 3-9 

327-12210-1- 

37.19* 

7. 5 2£ 

c 

1 9**171 

.. 3?7-l22i£-2 

.37.1.9* 

J*_ 52E 

c .. 

1 vk.ap? 0 

327-12220-1 

-3.25k 

6.03E 



32 7 - 4 2220-2 


6.C3E 


* ; ••• a > •* 9 

J5*»-1 226C-* 

. 7 : . 

*.17L 


* *) V A » 7 % 

327-1 2267-2 


. 1.'7E _ 


2 4 v A 1 - C 

’23-1 2790-* 

1 6 . 6 1 * 

3.2Fr 

c 

J 7 v A C 70 

3?f-1 2397-3 


3. ?8r 

c 

2C w A=!79 

22P-124O0-1 

* r . l ' k 

* .4PF 

c 

?*K?'5 

3? £-*2*7 0-2 

4 7 . 4 7t 

* . 4 C E 

r 

w 

??mA7'*9 

32 e -l 2*20-1 

- c . c * ' 

. 2 5 E 

c 

J 4m t c -9 

725-1 2*20-2 

*S.S*K 

.25* 

r 

2 1 V A “ 7 / 

3? 9-12570-1 . 

.561 

1 .43*' 


7 * ^ A e '7 c 

i?c-i 2570-2 

5 6 .5*» 

* .-3- 


?1 ^ AR-9 

329-12590-1 

45 . ClK 

3. 37w 


2.7 v t c 79 

t 29-1?5co-2 


3.37* 


21 vAR?9 

22 *>-i 3 no-i 

5* . 44K 

5 . 5 1 *• 


; 1 y A • 7 9 

3 29-1 3 0* 0-2 


5.51W 


23mA»79 

331-11 5P0-1 
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